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Abstract 

We demonstrate that the exact form-factors of the Calogero-Sutherland 
model which were recently found in in confirmation of the congectures 
earlier made by Haldane (Ref. Q) can be reproduced in the framework of 
some bosonization procedure in momentum space. This observation implies a 
possibility of an exact bosonization of this model describing one-dimensional 
anyons. 
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The Calogero-Sutherland (CS) model of one- dimensional (ID) spinless particles with 
long-ranged 1/x 2 interaction ( V(x) = ^ 2 p—wj^ i n a finite size L system) provides a remark- 
able example of an exactly solvable system which can be viewed as a continuous deformation 
away from either ideal Fermi or Bose gas. 

At special values A = 1/2, 1 and 2 the original continuous model describes statistics of 
eigenvalues of random matrices belonging to one of the three Dyson's ensembles. Recently 
these special cases were intensively studied in the framework of the theory of universal 
spectral correlations in random systems f2|. The lattice version of the CS model at A = 2 also 
induced a formulation of the Haldane-Shastry spin chain model [[J possessing rich Yangian 
symmetry in the case of a finite length chain Q . Further generalizations introducing particles 
with internal degrees of freedom in both continuous and lattice versions led to a family of 
new solvable models with graded SU(m\n) symmetries ||. 

The spectrum of the CS model obtainable by means of the so-called Asymptotic Bethe 
Ansatz (ABA) |J appears to be strikingly simple and enables one to consider this system as 
an ideal gas of excitations obeying fractional exclusion principle formulated in J7|. The latter 
was originally proposed as an alternative definition of fractional quantum statistics which 
is independent of the space dimension in contrast to the conventional notion of exchange 
statistics based on braiding relations for corresponding operators [§]. 

Remarkably, in the framework of the CS model these two definitions of fractional statis- 
tics demonstrate their perfect consistency. In its fractional exclusion formulation the prin- 
ciple of fractional statistics manifests itself as a set of rules governing allowed patterns of 
occupation numbers of states with momenta ki obeying ABA equations || for a system with 
periodic boundary conditions 

Lki = 2irrii + 7r(A — 1) sgn{ki — kj) (0.1) 

Eigenvalues of the Hamiltonian and the total momentum operator are given by free particle's 
expressions: E = ^Y,i and P = Y,iki- In the case of a rational coupling constant A = |C 
the ground state [Z{ = e~ ~) 
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*o({^}) = n(^-^) A iK A " 2 ~ (0.2) 

i<j i 

can be characterised as a periodic array of one occupied and M — 1 empty states associated 
with consequent non- negative integers rij. Patterns describing excited states are also sub- 
jected to a very restrictive selection rule which dictates that one can only move N momenta 
out of the ground state distribution simultaneously. This rule can be naturally interpreted in 
terms of fundamental quasiparticles with disperion e q {u) = ^(u 2 — kp), (\u\ > hp = ixp where 
p is a density of particles of mass m — 1) and quasiholes (eh( v ) = jx(kp — v 2 ), \v\ < kp) as 
a statement that an elementary excitation at conserved particle number is constituted by N 
quasiparticles and M quasiholes. On the basis of this constraint it was first anticipated in || 
and then explicitly found in [|K|] that dynamical correlation functions of local operators such 
as a fundamental field operator ty(x) and density p(x) have exclusively simple structure. It 
stems from the fact that sums over intermediate states \v > in formal expansions (where 
= %p) 

< vac\0(x, t)O(0, 0)\vac >= 

J2\ < v\O(0,0)\vac> | 2 exp(iE u t - iP v x) (0.3) 
now amount to single terms given by the only nonzero form-factors 



< ui, ...Miv_i;i;i, ...v M \ty(0)\vac >= CW-i 



nf (i - vf) 1 ^ nf- 1 ^ 2 - 1)- nJT" 1 ^ - < 



(0.4) 

and 



< ui, ...mjv; v u ...VM\p{0)\vac >= C m ,n —[- x — & -1=175 



nf K 2 - 1)^ nf (i - vf)^ nf M (n, - , 



(0.5) 

where C nm are known normalization factors |J, [|Kj and rapidities are normalized with 
respect to the Fermi momentum (A;, — > ki/kp). 

Formulae (4) for A = 2 [IT] and (5) for A = 2, 1 and 1/2 |T2|] were first obtained by means 



of supermatrix model calculus. In pi they were found at arbitrary rational A by using 



the method of Jack polynomials which was developed in [|T^] to compute static correlations 
functions related to Selberg integrals. 

This basically combinatoric approach explicitly implements the fractional exclusion prin- 
ciple. On the other hand, in view of the abovementioned equivalence of two definitions of 
fractional statistics one might expect that the problem also allows a purely algebraic solution 
based on monodromy properties of quasiparticle field operators and their braiding relations. 

Indeed, one might be tempted to interpret formulae (4)and (5) as vacuum expectation 
values of some "vertex operators" constructed from a free boson field ( 2D Coulomb gas 
correlators). In this Letter we present arguments in favor of such an interpretation which 
should be actually viewed as a possibility to carry out an exact bosonization of the CS 
model. 

In its original formulation the conventional bosonization procedure was developed to 
describe longwavelength properties of gapless ID fermions with contact interactions in terms 
of their collective excitations (charge and spin density waves). The procedure involves a 
coordinate space construction of the fundamental fermion field operator in terms of bosons 
(i/j(x) ~ Y,r,l e % ^ R ' L ^ ±%kFX ) which is suitable only for a description of low-energy fermionic 
excitations with linear dispersion in the vicinity of well separated Fermi points. Apparently, 
intending to develop an exact bosonization scheme valid at all scales one can no longer use 
this simple representation. 

Conceivably, some hints for a proper construction in coordinate space can be inferred 
from the fact that the wave function of the ground state (2) (as well as lowest excited states) 
can be reproduced as vacuum expectation values of a product of Af " quasi- chiral" vertex 
operators ip{z) = e 1 ^^) written in terms of some formal bosonic variable 

M 

*„({*}) =< jj e ^(**) e -<£CAr-i)*(o) > (0.6) 

i 

where the gaussian average has to be taken as < (f>(z)(f)(z') >= — \n(z — z'). Analogously, the 
wave function ^(z) = ^o(z) Ilf I\ J j' h (zi — Q) corresponding to the state with Afu quasiholes 
can be readily obtained by inserting a product of Afh operators e^^^ into (6). 



These observation are equivalent to those made about Laughlin trial wave functions f!4 
which are nothing but a result of analytic continuation of the CS model wave functions 
defined on a unit circle into an entire complex z plane (if nonanalytic gaussian factors 
e -N 2 / 4 ar6; as usual, asigned to integration measure). 

To make this heuristic construction physically meaningful we implement another obser- 
vation from the so-called collective field theory of unitary matrix models [|15|]. It was shown 



in 



15| that the system of free nonrelativistic fermions is equivalent to the cubic Hamiltonian 



h = |((^) 2 + (def) + \{d v fde + l -{d v f (0.7) 

where bosonic variables <p(x) and 9(x) are conjugated to each other: [<f(x), 0(y)) = 
Ysgn(x — y) and related to local particle number (p = ^dip) and current (j = ^09) den- 
sities. Elaborating on this theory one can see that although time evolution of bosonic 
variables is quite complicated, it turns out that some exponential composite operators 
W a b = exp(ia(p + ibO) may have much simpler dynamics. In particular, the exponentials 
Wi t x = exp(iif + i9) satisfies the equation of motion for free fermions. 

Moreover the exponentials Wa,i = exp(«A<p + i9) obey the equation of motion of the 
CS model at arbitrary A (modulo some normal ordering ambiguities). We believe that 
this circumstance makes it possible to consider W\,i as a plausible representation of the 
fundamental CS field with time dependence resulting from the evolution of variables (p(x) 
and 9{x) determined by the cubic Hamiltonian (7). It implies that the " quasi-chiral" boson 
field appearing in (6) can be identified as (f>(z) = Xip + 9. 

Referring to as a bare fermion ip(x) one can also treat as a "dressed" anyon 
field operator: 

$(x) = e^-^r^^v^v^) (0.8) 

which is another formulation of the Eq.(l). In terms of ip(x) the CS Hamiltonian becomes 
quadratic which sheds additional light on the remarkable simplicity of the spectrum of the 
CS model describing purely statistical interaction. 



However this real space bosonization approach doesn't provide an immediate construction 
of eigenstates of the model. Even for the case A = 1 studied in the context of the 1 + 1- 



dimensional string theory in [16j it appears to be quite involved. 



On the other hand, this aim can be readily achieved in the framework of a dual ap- 



proach of "momentum space bosonization" recently proposed in [I7|. Being applied to the 
Sine-Gordon model at (5 — 1 this method provided a bosonic representation of a massive 
relativistic fermion field which facilitates computations of formfactors of other operators 
built from the Sine-Gordon field. This study also clarified a physical origin of the bosonic 
field <p(z) from (6) as a generating function for the integrals of motion I n = YliPi of the free 
fermion system 

4>{z) = £ - (0-9) 
"I in 

where I n at n < are "angle" variables conjugated to those at n > ("action" variables). 
Altogether they form a U(l) Kac-Moody algebra: [I n , I m ] = nS nm . 

Being formulated this way the method of |T7| allows a straightforward generalization 
onto the case of interacting systems which are characterised by factorizable scattering and 
possess M integrals of motion given by powers of some Lax operator: I n = TrL n . Since 
the CS model determined by the Lax operator = 5ijPi + ^(1 — Sij)cot t{x% — Xj) allows 
the ABA solution, the eigenvalues of integrals of motion can be also presented in the free 
particle form I n \{ki} >= J2ik^\{ki} > where fc, are asymptotic momenta obtainable from 
(6). 

In the picture of second quantization the integrals of motion acquire free particle form 
in terms of quasiparticle and quasihole creation- annihilation operators Z a (u) (a = p, h) 
labeled by asymptotic momenta. Being considered on general grounds operators Z a obey 
Zamolodchikov-Faddeev (ZF) algebra (a = 1, ...,/) |p~8| : 

Z a (u)Z b (v) = S ab (u - v)Z b (v)Z a (u) 
Z\{u)Z\{v) = S ab (u-v)Zl(v)Zl(u) 
Z a {u)Z\{v) = S~ b \u - v)zl(v)Z a (u) + 6 ab 5(u - v) (0.10) 



where S a b(u — v) is a scattering matrix of excitations which depends on relative rapidity. In 
the case S ab (u) = ^albSiu) is a U(l) "charge" of a species a) bosonization procedure 
can be performed by constructing a representation of the algebra (10) Z a (u) = e na$<M ) in 
terms of a one-component bosonic field with commutation relations 

[$(«),$(«)] =]nS(u-v) (0.11) 

which infer the vacuum expectation value 

< 0|$(u)$(u)|0 >= -lnf(u-v) (0.12) 

with f(u) determined by the equation S(u) = f(—u)/f(u). In terms of ZF operators the 
integrals of motion can be written as 

/ , 

/„ = E / duZl(u)p n a (u)Z a (u) (0.13) 

where A a is a region of rapidities supporting excitations Z a with momenta p a (u). 

Expressions for local (in momentum space) operators bilinear in Z a can be found in a 
manner similar to the point-splitting procedure in the conventional real space bosonization 
which yields 

*mz.m = f-^> (o.i4) 

With the use of (14) the bosonic phase (9) providing a real space representation of the 
fundamental field operator ip(z) can be recovered in the form 

0W = E-/ — ln(l-«V)^ (0.15) 
vy ! V 2 Jr a 2ni K J du K J 

where integrals in the complex plane of rapidities are taken over contours T a encircling real 
axis intervals A a . 

In the particular case of the CS model scattering matrix is given by matrix elements 
Spp(ui - u 2 ) = ir\sgn(ui - u 2 ), S hh (vi - v 2 ) = \sgn{vi - v 2 ) and S ph (u - v) = 7isgn(u - v) 
equal to statistical phase factors assigned to elementary excitations having 7^ = iA^ 1 / 2 
with rapidities equal to their velocities. It then leads to ZF operators 



Z p {u) = e^ u \ Z h (v) = e- i *% (0.16) 

creating eigenstates (\uf, Vj >= flj ^p(*0 Uj ^l( v j)\ vac >) an d < 0|$(-u)$(t>)|0 >= — ln(-u — 
f). By construction is subjected to the condition $(0) = and has no singularities 
on the real axis. Naturally, $(«) should also vanish at \u\ — > oo. Although general rules 
of choosing integration contours in (15) remain to be clarified we believe that a correct 
prescription is to take contour intergals at positive and negative Rew separately to keep a 
distinction between left and right moving excitations. 

Then integrals taken over four contours encircling intervals (— oo, —1], [—1, 0), (0, 1], and 
[1, oo) yield the following expression for the fundamental CS field operator at z — 1 (x — 0): 

^(0) ~ exp^/A- -^=)($(1) + $(-1)) (0.17) 

Together with the formulae for Z p (u) and Zh(v) given by (16) the gaussian average 

N-l M 

<u 1 ,...u N - 1 ;v 1 ,...v M \'^(0)\vac>=< JJ Z p iui)J\_Z h {vj)^) > (0.18) 

i j 

reproduces the form- factor (4) up to a normalization constant. The abovementioned se- 
lection rules constraining allowed quasiparticle contents of asymptotic states provide that 
(18) is free of infrared divergent factors potentially coming out of averages < <3> 2 (-u) > at 
coincident momenta. 

Notice that one can only take the free fermion limit (M = N — 1) in the final formula 
(4) to obtain < v\ip(0)\vac >= 6(1 - \v\). 

By inspection it is also easy to see that to reproduce the form-factor of the density 
operator (5) one can use the following prescription 

p(0) ~ Pexp^VX- -^=)($(1) + $(-1)) (0.19) 

where P = £fui-j£f Vj is a total momentum of excitations which can be identified with 
(d z (f)(z) — dz4>(z))\ z ^i. The relation (19) should result from an operator algebra of ip(z) and 
fjy{z) at arbitrary z which remains to be established. 
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In conclusion we notice that a systematic approach to a complete bosonization of the 



CS model could be built on the basis of its infinite symmetry recognised as [yj . The 
corresponding linear algebra is formed by operators W£ = Yn' xl~ x p^ +a ~ l + ... which can be 
recursively obtained as = 2 ^ x +s ^ [Yjf x h Wn+2] from the CS integrals of motion: W x = I n . 
In particular, the CS Hamiltonian is among these operators (Hcs = W^)- Representation 



theory of in terms of free bosons discussed in [^C§ might appear to be useful here. 

A complete CS model bosonization scheme valid at all scales would be an important 
step towards an exact formulation of a (nonlinear) hydrodynamics of ID fluids of interacting 
fermions. 

The author is indebted to F.D.M.Haldane for valuable discussions and for the opportunity 



to become aware of Refs. || and |L0| prior to publication. This work was supported by the 
NSF Grant. 
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